■’^’Calhoun 

IflitiKuliarttl  Arehivtaf  (he  Naval  ftutgrdduH*  School 


Authors) 

Nicholson,  Harry  E. 

Title 

An  analysis  of  a  divergent  model  for  numerical  forecast]  ng. 

Publisher 

Monterey,  California.  Naval  Postgraduate  School 

Issue  Date 

1962 

URL 

http://hdl  .handle  ndt/10945/11607 

This  docunnent  was  downloaded  on  April  24,  2015  at  04:54:20 


DUDLEY 

KNOX 

LIBRARY 


htt  p:/ /www.  n  ps.  e  d  u/l  ib  ra  ry 


Calhoun  is  a  project  of  the  Dudley  Knox  library  at  NP5,  furthering  the  precepts  and 
goals  of  open  government:  and  government  transparency.  All  information  contained 
herein  has  been  approved  for  release  by  the  NPS  Public  Affairs  Officer. 

Dudley  Knox  Library  /  Naval  Postgraduate  School 
411  Dyer  Road  /  1  University  Circle 
Monterey,  California  USA  93943 


http://www.n  ps .  edu/ 


v  ■:•  ..  .>r. •’i  ••“"  "•  •"  '■.%%  '-V-  vk * 

>«  .■',  JLc,-*  -'  \.  -•■-  "  ■..  ’•£  I  ■•'  ,  •  •«.  *- 

•::.V  .  /  «-v^;vmc: 


AN  ANALYSIS  OF  A  DIVERGENT  MODEL 
'  FOR  NUMERICAL  FORECASTING 


HARRY  E.  NICHOLSON 


NPS  ARCHIVE 
1962 

NICHOLSON,  H. 


AN  ANALYSIS  OF  A  DIVERGENT 
MODEL  FOR  NUMERICAL  FORECASTING 


-*  # 


Harry  E.  Nicholson 


AN  ANALYSIS  OF  A  DIVERGENT 


MODEL  FOR  NUMERICAL  FORECASTING 


by 

Harry  E.  Nicholson 

!t 

Lieutenant,  United  States  Navy- 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 

METEOROLOGY 


United  States  Naval  Postgraduate  School 
Monterey,  California 


19  6  2 


library  DUDLEY  KNOX  LIBRARY 

U  S.  NAVAL  POSTGRADUATE  SCHOOL  NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA  MONTEREY,  CA  93943-5101 


AN  ANALYSIS  OF  A  DIVERGENT 
MODEL  FOR  NUMERICAL  FORECASTING 


by 

Harry  E.  Nicholson 


This  vjork  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 

METEOROLOGY 
.  from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 


A  forecast  model  proposed  by  'Arnason  as  being  capable  of  producing 
realistic  forecasts  of  displacement  of  systems  at  and  below  500  mb  is 
investigated.  The  prediction  equation, 


is  applied  to  test  cases  at  1000,  850,  700  and  500  mb  with  various  values 
of  the  parameter  k]_.  The  effect  of  the  value  of  the  parameter  on  the 
amount  of  divergence  present  in  the  model  is  discussed.  Root -mean-square 
errors  for  2lj-hour  and  ij8-hour  forecasts  with  various  values  of  k^  at 
each  level  are  presented.  Optimum  values  of  k-^  for  the  cases  tested  are 
selected  as  3.5*  2.?5  and  0.9  for  the  levels  850,  700  and  500  mb  respec¬ 
tively.  It  is  found  that  with  proper  values  for  k-^  the  model  will  provide 
good  forecasts  of  system  movement  at  each  level  tested,  but  as  expected, 
changes  in  the  intensity  of  systems  are  not  correctly  forecast.  In  all 
cases  values  of  root-mean- square  error  of  forecasts  are  less  than  those 
obtained  by  applying  the  "Helmholtz"  barotropic  model.  A  geostrophic 
version  of  the  model  is  also  tested  and  is  found  to  produce  results  simi¬ 
lar  to  those  obtained  for  the  stream-function  version  except  that  the 
geostrophic  "blow-up"  of  low-latitude  high  pressure  cells  is  observed. 

The  writer  wishes  to  express  his  appreciation  for  the  guidance  and 
encouragement  given  him  by  Professor  George  J.  Haltiner  of  the  U.  S. 
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tion,  appreciation  is  expressed  to  Mr.  Geirmundur  ’Arnason,  on  whose  work 
this  paper  is  based,  and  to  personnel  of  the  Fleet  Numerical  Weather 
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Introduction 


Since  the  initation  of  operational  numerical  weather  forecasts  by 
the  Joint  Numerical  Weather  Prediction  Unit  in  1955.,  the  only  forecast 
model  which  has  stood  the  test  of  daily  operational  use  has  been  the 
barotropic  model  first  proposed  by  Rossby  in  1939.  The  barotropic  model 
remained  in  its  original  form  until  1958  when  Wolff  [1]  discovered  large 
systematic  errors  due  to  improper  forecast  of  very  long  waves.  After 
these  errors  had  been  empirically  corrected  by  Wolff,  it  was  discovered 
by  Cressman  [2]  that  these  long-wave  errors  could  be  reduced  by  the  in¬ 
troduction  of  divergence  into  the  barotropic  model.  This  resulted  in 
the  divergent  or  "Helmholtz"  version  of  the  barotropic  model  which  is  in 
general  use  at  this  time. 

With  the  improved  barotropic  model,  high  quality  forecasts  have  been 
available  for  the  500-mb  level.  This  barotropic  model,  however,  is  by 
nature  restricted  in  its  use  to  the  500-mb  level  and,  in  addition,  is  ca¬ 
pable  mainly  of  forecasting  the  movement  of  systems,  and  not  in  general 
changes  in  intensity.  Occasionally  attempts  have  been  made  to  replace 
the  barotropic  model  in  operational  use  with  some  type  of  baroclinic 
model  which  would  obviate  these  shortcomings.  These  attempts  have  so  far 
met  with  failure  largely  due  to  the  inability  to  control  the  forecast  of 
cyclogenesis. 

’Arnason  has  suggested  that,  in  view  of  the  great  utility  of  the 
500-mb  forecasts  produced  by  the  barotropic  model,  it  might  be  worth¬ 
while  temporarily  to  set  aside  the  problem  of  predicting  cyclogenesis 
and  to  try  to  develop  a  model  which  would  possess  the  same  skill  of  fore¬ 
casting  movement  at  other  levels  which  the  barotropic  has  at  500  mb.  To 
this  end,  ’Arnason  has  proposed  a  model  [ 3 i]  which  maintains  many  of  the 
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characteristics  of  the  barotropic,  but  which  is  more  adaptable  to  other 
levels. 

This  model  is  based  on  a  single-layer  compressible  fluid  with  a 
stratified  density  distribution,  bounded  by  two  constant-pressure  sur¬ 
faces.  The  theoretical  development  results  in  an  expression  for  the 
horizontal  wind  divergence 

JwW*  j +  V-Vs»\  (i) 

J'HUi  /  . 

Introducing  this  expression  into  the  vorticity  equation, 

+  V- V7-7  d>v  V  s  O  j  (2) 

gives  the  prediction  equation  for  the  "stratified"  model 

/If +V'7A  =  °.  (3) 

.H  '  / 

In  order  to  apply  this  equation  to  forecasting,  'Arnason  has  suggested, 
as  one  alternative,  the  use  of  a  stream  function  for  computation  of  wind 
and  vorticity  and  the  elimination  of  <f>  through  a  geostrophic  assump¬ 
tion.  This  gives 

/ ^  11  -v  w-v-y  -  r 7  v  •  7  ?  -  0 ,  a) 

v  ?H/d*  /  *'H 

which  is  of  the  same  form  as  the  prediction  equation  for  the  divergent 
or  "Helmholtz"  barotropic  model,  but  with  an  additional  divergence  term 
involving  V  and  with  the  substitution  of  the  stability-dependent  term 
g*  for  g. 

Although  the  forms  of  the  prediction  equations  are  similar,  there 
are  several  basic  differences  between  the  "stratified"  and  other  barotrop¬ 
ic  models.  By  the  method  of  perturbation  analysis  'Arnason  has  obtained 
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the  meteorologically  significant  root  of  the  frequency  equation  for  the 
stratified  model 


C  = 


(5) 


Although  this  equation  is  again  similar  to  the  corresponding  equation 
for  the  Helmholtz  barotropic  model,  an  important  difference  arises  in 
the  fact  that  the  basic  current  u  does  not  have  a  multiplying  factor. 
By  empirically  varying  the  value  of  g'H,  the  magnitude  of  the  ^  term 
may  be  changed  without  reducing  the  effect  of  the  basic  current,  thus 
allowing  control  over  the  speed  of  movement  of  systems.  As  can  be  seen 
from  equation  (l)  this  is  equivalent  to  a  variation  in  the  amount  of 
divergence  present  in  the  model.  This  provides  a  method  of  adjusting 
the  model  to  yield  proper  system  movement  at  a  number  of  levels. 

The  stratified  model  as  proposed  by  'Arnason  seems  to  provide  a 
promising  approach  to  numerical  forecasting  at  and  below  the  500-mb 
level.  It  is  the  purpose  of  this  paper  to  apply  the  model  to  a  number 
of  test  cases  and  determine  its  effectiveness. 
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Application  of  the  Model  to  Forecasting 

In  the  application  of  the  prediction  equation  (3)  to  actual  fore¬ 
casts  there  are  several  choices  to  be  made.  The  parameter  V  has  not 
been  explicitly  defined,  but  may  be  taken  to  be  any  kind  of  space- 
averaged  or  zonal  wind.  Also  the  value  of  gfH  is  left  to  be  determined 
empirically. 

By  examination  of  the  prediction  equation  (3)*  it  can  be  seen 
that  the  additional  ''divergence"  term  of  the  stratified  model  which  in¬ 
volves  V  essentially  represents  the  horizontal  advection  of  height 
by  the  mean  wind.  This  immediately  suggests  the  notion  of  a  steering 
current,  which  may  be  thought  of  as  a  wind  at  some  level  or  as  a  mean 
wind  in  some  layer  of  the  atmosphere.  Also,  the  magnitude  of  the  di¬ 
vergence  term  is  dependent  on  the  difference  between  the  mean  wind  and 
the  wind  at  the  forecast  level.  Thus  it  appears  reasonable  to  use  as 
a  mean  wind  something  close  to  the  5>00-mb  wind,  since  this  would  mini¬ 
mize  the  divergence  term  at  500  mb  and  allow  it  to  increase  with  dis¬ 
tance  from  500  mb.  This  is  generally  consistent  both  with  observed 
patterns  of  divergence  in  the  atmosphere  and  with  the  need  for  increas¬ 
ing  the  value  of  the  divergence  term  at  lower  levels  to  supplement  the 
movement  due  to  vorticity  advection  alone. 

A  second  parameter  which  must  be  empirically  interpreted  is  the 
factor  g'H.  The  value  of  g1  is  a  measure  of  the  static  stability  of 
the  basic  flow  and  H  is  the  thickness  of  the  layer  considered  in  the 
model.  Since  neither  is  uniquely  defined,  the  combination  may  be  taken 
as  any  reasonable  value  which  produces  proper  movement  of  the  systems 
being  forecast. 

In  order  to  make  a  forecast,  either  V  or  $  must  be  eliminated 

h 


from  equation  (3).  'Arnason  has  used  a  geostrophic  relation  to  eliminate 
^ ,  obtaining  equation  (I4 )  which  involves  f  alone.  To  apply  this  equa¬ 
tion,  height  fields  must  first  be  "balanced"  to  obtain  the  stream  function 
and  then,  at  the  conclusion  of  the  forecast,  inverted  to  obtain  the  fore¬ 
cast  height  field.  Another  approach  would  be  to  compute  the  vorticity 
and  winds  directly  from  the  height  values  using  the  geostrophic  rela¬ 
tionship,  obtaining  the  forecast  equation; 

(VK£X\  ^  t  p. \V'V7-£l  V  ■  V2  =  o  (6) 

v  pT/IT  %  '  fH 
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3.  Procedure 


As  mentioned  above,  in  order  to  apply  the  stratified  model  to  actual 
forecasting,  values  must  be  chosen  for  the  parameters  V  and  g'H.  For 
the  purposes  of  the  tests  presented  here,  it  was  decided  to  use  for  \/ 
a  mean  wind  computed  from  a  height  field  obtained  by  first  heavily 
smoothing  the  500-mb  and  700-mb  heights,  then  computing  an  intermediate 
field  by  the  equation 

^  =  >1S  +  >2  5  $  (7) 

500  *700 

The  wind  obtained  by  this  method  is  thus  a  space-mean  wind  for  some 
intermediate  level  between  500  mb  and  700  mb.  For  the  first  six  one- 
hour  time  steps,  a  V  computed  from  the  initial  data  was  used;  there¬ 
after  a  new  ^7  was  computed  from  forecast  fields  every  sixth  time  step. 
In  the  choice  of  the  parameter  g'H,  consideration  was  given  to  the 
similarity  between  the  'forecast  equation  (U)  and  the  forecast  equation 
for  the  "Helmholtz"  barotropic 


Since  it  was  known  that  forecasts  made  with  this  model  at  500  mb  are  of 
consistently  good  quality,  it  was  decided  to  choose  g'H  such  that  the 
Helmholtz  term  for  the  stratified  model  would  be  the  same  as  the  Helm¬ 
holtz  term  of  the  barotropic.  Thus  in  the  Helmholtz  term,  g'H  was 
chosen  to  be 

'H  *  o  M.  (9) 

where  k  is  a  "tuning  constant"  with  a  value  of  four,  and  H^  is  the  stand¬ 
ard  height  of  the  500-mb  surface. 

In  order  that  the  model  be  used  at  levels  other  than  500  mb,  it  was 
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decided  to  include  a  second  "tuning  constant"  in  the  final  divergence 
term.  By  varying  this  second  constant  the  proper  amount  of  divergence 
may  be  introduced  at  each  level  to  produce  the  best  forecast.  With 
these  values  introduced,  the  forecast  equation  for  the  stream  function 
case  becomes 

(  7  V'Vt-"0  (10) 

In  the  same  manner,  starting  with  equation  (6)  a  similar  equation  was 
obtained  for  the  geostrophic  version  of  the  model: 

/ v*-±£iA  Si  +£V'V7 - iuJt.  v- V2:  *  o  (n) 
'  gHT/  JT  %  ( 

The  model,  in  the  form  of  finite  difference  equations  obtained  from 
equations  (10)  and  (ll),  was  run  to  lj8~hour  forecasts  for  850  mb,  700  mb 
and  500  mb  for  several  randomly  selected  days.  Various  values  of  k-^  were 
used  at  each  level  for  both  the  geostrophic  and  stream-function  cases. 

In  addition,  for  the  6  May  1962  case  runs  were  made  on  1000-mb  data  in 
order  to  determine  the  feasibility  of  using  the  model  at  that  level. 
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h.  Computational  Methods 

The  test  computations  of  the  stratified  model  described  in  this 
paper  were  made  on  a  Control  Data  Corporation  l60l|.  computer  using  a 
program  constructed  of  meteorological  subroutines  prepared  by  the  Fleet 
Numerical  Weather  Facility,  Monterey.  The  grid  used  in  the  numerical 
solution  of  the  prediction  equation  was  the  standard  octagonal  grid  used 
by  both  the  Joint  Numerical  Weather  Prediction  Unit  at  Suitland,  Md.  and 
the  Fleet  Numerical  Weather  Facility. 

Time  differencing  was  accomplished  in  one-hour  time  steps  using  the 
central  difference  formula 

i  +  2(M)y  ,  (12) 

where  A  +  4*  is  the  computed  one-hour  change  in  geopotential.  At  the 
initial  time,  T-0  ,  and  all  integral  multiples  of  twelve,  the  forward 

difference  formula 

=  K  +  (4*  Dr  (13> 

was  used. 

In  cases  where  a  stream  function  wind  was  used,  the  stream  values 
were  obtained  from  initial  height  data  through  the  solution  of  the  bal¬ 
ance  equation  as  described  by  'Arnason  All  initial  height  fields 

and  fields  required  for  forecast  verification  were  obtained  from  files 
of  operational  analyses  at  the  Fleet  Numerical  Weather  Facility. 
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5.  Results 

Verification  scores  for  each  of  the  cases  tested  were  determined 
by  a  FNWF  program  using  the  equations 


aw 


Pillow 


N 


(15) 


Here  A  represents  the  forecast  value,  B,  the  observed  value  and  N,  the 
total  number  of  grid  points,  1977. 

As  can  be  seen  from  the  equations,  the  pillow  is  a  measure  of  the 
error  in  the  forecast  of  the  mean  height  while  the  RMSE  (root  mean  square 
error)  is  a  measure  of  the  deviation  of  errors  around  the  mean  error. 

Figures  1,  2  and  3  show  plots  of  RMSE  as  a  function  of  k-^  for  the 
various  cases  tested.  As  can  be  seen,  the  value  of  k-^  necessary  to  pro¬ 
duce  the  best  forecast  (minimum  RMSE)  varies  somewhat  from  case  to  case 
at  each  level.  However,  the  curves  are  relatively  flat  in  the  vicinity 
of  the  minima  so  that  by  selecting  a  reasonable  value  of  k-p  the  quali¬ 
ty  of  forecasts  from  one  case  to  another  will  vary  only  slightly  from 
the  optimum.  Table  1  shows  mean  values  of  RMSE  for  all  cases  obtained 
with  various  values  of  k^  at  each  level.  A  plot  of  these  values  is 
shown  in  figure  i*.  From  this  plot  of  mean  error  versus  k-p  values  of 
kj_  can  be  selected  which  will  give  the  minimum  total  error.  Table  2 
gives  a  tabulation  of  the  error  which  would  occur  in  each  case  tested  if 
kj  at  each  level  were  selected  as  the  value  which  produces  the  smallest 
mean  error. 

Forecasts  for  850  mb  made  from  initial  data  for  20  January  1962  are 
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Table  1.  Mean  Verification  Scores  of  Forecasts  With  Various  Values  of 


STREAM-FUNCTION  FORECASTS 

ki 

RMSE  (ft) 

kl 

RMSE  (ft) 
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RMSE  (ft) 

CO 

850 

2 

112 

3 

107 

h 

108 

5 

lilt 

u 

t-4 

p 

o 

700 

1 

133 

2 

12U 

3 

123 

h 

128 

pd 

C\J 

5oo 

0 

12*8 

1 

2 

1U7 

1 

1U7 

2 

150 
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SI 

700 

1 
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2 
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h 
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1 
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1 
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2 
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j 
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GE0STR0PHIC  FORECASTS 

t 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

) 

CO 
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2 
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3 
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h 
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5 
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U 

Jd 

o 

s: 
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2 
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♦ 

3 
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h 

131 
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0 
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1 

2 
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2 
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CO 

850 

2 

198 

3 

192 

h 

193 

5 

200 

u 

!d 

o 
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700 

1 

— 

2 

203 

3 

195 

k 

198 

oo 

pd- 

5oo 

0 

1 

2 

271 

1 

269 

2 

277 
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Table  2.  Verification  Scores  For  Each  Day  Using  Optimum  Value  of 


850  mb 

RMSE  (ft) 

700  mb 

RMSE  (ft) 

500  mb 

RMSE  (ft) 

k]_=3.5 

kx=2.75 

k^O.9 

1*  November  1961 
(2k  hr) 

111 

116 

1 31* 

20  January  1962 
(21*  hr) 

128 

128 

178 

6  May  1962 
(2k  hr) 

87 

88 

127 

1*  November  1961 
(1*8  hr) 

173 

185 

226 

20  January  1962 
(1*8  hr) 

198 

185 

250 

6  May  1962 
(1*8  hr) 

11*3 

l!*0 

168 
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shown  in  figures  8  through  13.  Initial  charts  and  charts  for  the  veri¬ 
fying  time  of  each  forecast  are  shown  in  figures  5  through  7.  Comparing 
forecasts  for  the  various  values  of  k-^  at  a  given  level  shows  that  the 
forecast  speed  of  synoptic  waves  increases  as  the  value  of  is  increased 

(thus  decreasing  the  value  of  g’H  and  increasing  the  magnitude  of  diver¬ 
gence).  This  is  in  accordance  with  the  results  of  equation  (5).  Further, 
it  can  be  seen  that  there  is  very  little  change  in  intensity  of  the  vari¬ 
ous  systems  as  k-^  is  varied.  This  demonstrates  the  point  that  the  addi¬ 
tional  divergence  term  of  the  stratified  model  results  in  movement  rather 
than  development. 

As  can  be  seen  from  equations  (10)  and  (ll)  a  choice  of  k^=0  would 
reduce  the  forecast  to  the  barotropic.  An  examination  of  figure  re¬ 
veals  that  the  optimum  value  of  k]_  at  500  mb  is  close  to  zero,  but  that 
the  stratified  model  does  show  a  slight  improvement  over  the  barotropic 
with  values  of  k^  between  zero  and  one.  This  is  as  expected,  since  500 
mb  is  generally  accepted  to  be  at  or  near  the  level  of  minimum  divergence, 
and  thus  the  contribution  of  divergence  to  the  change  in  the  vorticity 
pattern  at  this  level  should  be  small. 

At  levels  other  than  500  mb,  where  the  barotropic  model  is  less  ade¬ 
quate,  the  values  of  k-j_  for  which  the  minimum  RMSE  is  reached  are  greater. 
For  the  700-mb  level  optimum  values  for  the  cases  tested  range  from  k]®2 
to  ki=Uj  while  at  850  mb  they  range  from  3  to  5.  These  results  can  be 
accounted  for  by  the  vertical  distribution  of  velocity  divergence  which 
generally  increases  with  distance  from  the  500-mb  surface. 

By  use  of  equation  (1)  and  the  chosen  values  for  g’H  given  by  equa¬ 
tion  (9)*  an  expression  for  the  divergence  at  a  given  time  may  be  ob¬ 
tained: 
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(16) 


cKvV  -  q-P  [11  +  k,\v- Vf 

In  figures  llj,  13>  and  16  divergence  charts  at  8£0,  700  and  $00  mb  as  com¬ 
puted  from  equation  (16)  for  initial  data  of  20  January  1962  are  shown. 
These  charts  show  clearly  the  vertical  distribution  of  divergence,  with 
values  as  large  as  +  5  x  10"^  sec-"*"  at  85>0  mb,  +  3  x  10 sec  ^  at  700 
mb  and  values  very  nearly  zero  at  f?00  mb. 

An  examination  of  the  divergence  fields  relative  to  their  correspond¬ 
ing  height  fields  reveals  that  lines  of  zero  divergence  very  closely  fol¬ 
low  the  short-wave  trough  and  ridge  lines  of  the  corresponding  height 
fields.  Between  the  short-wave  trough  and  ridge  lines  are  located  the 
centers  of  maximum  convergence  and  divergence,  with  divergence  ahead  of 
ridges  and  convergence  ahead  of  troughs,  as  might  be  expected. 

For  the  case  of  6  May  1962,  1000-mb  forecasts  were  made.  From 
figure  3  it  can  be  seen  that  the  curves  of  RMSE  versus  k]_  for  1000  mb 
are  almost  exactly  parallel  to  the  curves  for  8£0  mb  with  an  RMSE  about 
ten  feet  greater  at  1000  mb  than  at  85>0  mb.  The  fact  that  the  curves 
are  nearly  parallel  suggests  that  the  increased  error  at  1000  mb  is  due 
to  a  greater  change  in  the  intensity  of  systems  at  this  level  than  at 
8£0  mb.  The  movement  of  established  systems  which  appear  on  both  initial 
charts  is  nearly  the  same  at  both  levels ;  thus  k-^  which  controls  the  mag¬ 
nitude  of  the  additional  movement  term  should  reasonably  be  chosen  the 
same  for  the  two  levels.  The  verifying  charts  for  the  forecasts  in 
question,  figures  21,  22,  2E>  and  26,  show  this  to  be  the  case.  Forecasts 
at  both  levels  have  the  major  systems  placed  very  close  to  their  actual 
positions  at  verifying  time,  but  the  verifying  charts  at  1000  mb  differ 
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more  in  pattern  due  to  the  variation  in  intensity  change  of  different 
systems  over  the  forecast  period. 

In  addition  to  the  forecasts  made  for  the  three  days  using  stream 
functions  obtained  from  the  initial  height  fields,  for  two  of  the  days 
forecasts  were  made  directly  from  the  height  fields  using  equation  (11). 
As  can  be  seen  from  figures  27  through  30,  the  geostrophic  forecasts 
thus  obtained  showed  great  similarity  to  the  forecasts  made  from  the 
stream  function.  However,  as  expected,  the  geostrophic  version  appears 
to  suffer  from  the  major  error  common  to  geostrophic  forecasts  -  the 
"blow-up"  of  low- latitude  high-pressure  areas.  As  can  be  seen,  the  lj8- 
hour  geostrophic  forecast  increases  the  central  height  of  the  500-mb 
high  off  the  east  coast  of  the  United  States  to  over  19,500  feet,  while  “ 
the  stream-function  version  keeps  the  central  height  to  a  more  reasonable 
value.  The  treatment  of  the  high  off  the  west  coast  of  the  United  States 
is  similar. 
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6.  Conclusions 

The  stratified  model  as  proposed  by  ’Arnason  appears,  on  the  basis 
of  the  test  cases  presented  here,  to  produce  reasonable  forecasts  for 
levels  at  and  beloitf  500  mb.  The  stream-function  version  of  the  strati¬ 
fied  model  produces  500-mb  forecasts  for  2k  and  I48  hours  which  are  in 
all  cases  very  slightly  better  than  those  produced  by  the  barotropic 
model.  At  850  mb  and  700  mb  the  model  produces  forecasts  of  system 
movement  which  are  very  realistic.  The  geostrophic  version  of  the  model 
suffers  from  the  "blow-up"  of  low-latitude  high  pressure  cells,  as  is 
common  with  most  geostrophic  models. 

Optimum  values  of  the  multiplying  factor  of  the  divergence  term 
k]_  at  850,  700  and  500  mb  have  been  found  to  be  3.5,  2.75  and  0.9 
respectively.  Although  the  number  of  cases  tested  was  not  large,  the 
values  determined  here  should  be  sufficiently  accurate  to  provide  a 
basis  for  the  long-term  testing  necessary  before  any  model  becomes  suit¬ 
able  for  operational  use. 
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Figure  1;.  Plot  of  vs.  mean  RMSE 
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Figure  850-mb  analysis. 


20  January  1962 
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Figure  6.  8^0-rib  analysis,  21  January  1962 
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Figure  7«  8£0-mb  analysis,  22  January  1962 
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Figure  8. 


2U-hour  forecast  for  8j?0  mb  from  20  January  1962,  k-^0 


23 


Figure  9«  2l;-hour  forecast  for  8^0  mb  from  20  Januaiy  1962,  k.^-3 
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Figure  10.  2l*-hour  forecast  for  850  mb  from  20  January  1962,  k^=5 
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Figure  11.  l4.8-h.our  forecast  for  3^0  mb  from  20  January  1962,  k-^=0 
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Figure  12.  i|8-hour  forecast  for  850  mb  from  20  January  1962,  k^=3 
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Figure  13.  US-hour  forecast  for  8j?0  mb  from  20  January  1962,  k^=5 
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Figure  lU.  850-mb  divergence,  20  January  1962 
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Figure  1$,  700-rab  divergence,  20  January  1962 


30 


contour  interval 


2  x 


-1 

sec 


Figure  16.  j?00-mb  divergence,  20  January  1962 
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Figure  17.  1000-mb  analysis,  6  Kay  1 962 
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Figure  18.  850-iab  analysis,  6  May  1962 
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Figure  19.  2l*-hour  forecast  for  1000  mb  from  6  May  1962 
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Figure  20.  2li-hour  forecast  for  850  mb  from  6  Fay  1962 
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Figure  21.  1000-mb  analysis,  7  Kay  1962 


Figure  22.  850-mb  analysis,  7  May  1962 
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Figure  23.  U8-hour  forecast  for  1000  mb  from  6  May  1962 
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Figure  2U.  U8-hour  forecast  for  850  mb  from  6  May  1962 
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Figure  2 5*  1000-mb  analysis,  8  May  1962 
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Figure  26.  850-mb  analysis,  8  May  1962 
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b  November  1961 
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Figure  27 •  500-rob  analysis. 


Figure  28.  U8-hour  forecast  for  500  mb  from  U  November  1961 
(stream  function) 
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Figure  29 •  l|3-hour  forecast  for  500  mb  from  U  November  1961 
(geostrophic) 
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Figure  30.  500-mb  analysis, 


6  November  1961 
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